Origin of the T c enhancement in heterostructure cuprate superconductors 
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Recent experiments on heterostructures composed of two or more films of cuprate superconductors 
of different oxygen doping levels have shown a remarkable T c enhancement (up to 50%) relative 
to single compound films. We provide here a simple explanation of the enhancement which arises 
naturally from a collection of experimental works. We show that the enhancement could be caused 
by a structural change in the lattice, namely an increase in the distance of the apical oxygen from 
the copper-oxygen plane. This increase modifies the effective off-site interaction in the plane which 
in turn enhances the d-wave superconductivity order parameter. To illustrate this point we study 
the extended Hubbard model using the fluctuation exchange approximation. 



In 2008 Yuli at ai. [H presented a novel technique to en- 
hance the superconducting transition temperature (T c ) 
of LSCO films. The enhancement was achieved in a het- 
erostructure of two thin layers of the same cuprate super- 
conductor with different carrier density in the two layers. 
The experimental setup used I^-^Sr^CuC^ (LSCO) 
structures where the top layer had x — 0.35 (overdoped) 
and in the bottom layer x ranged from the underdoped 
to the overdoped. In all measured heterostructures T c 
was higher in the heterostructure than that of a single 
compound LSCO film with the same doping x as the 
bottom layer. The largest enhancement occurred around 
x = 0.12 where T c of the heterostructure was 32K, about 
50% larger than T c of the single compound film. 

In another experiment Gozar et aZ.[2j created het- 
erostructures with high T c from layers of metallic LSCO 
(overdoped) and either insulating LCO (the parent com- 
pound) or superconducting LCO (in this experiment the 
compound was LaCu04+a with 6 = for the insulator). 

These intriguing results present a possible new direc- 
tion in the exploration of high temperature superconduc- 
tors. At the moment, the phenomenon of T c enhance- 
ment is far from understood and seems to depend cru- 
cially on materials and growth method^. A few ideas 
have been put forth regarding the mechanism by which 
the enhancement occurs. These ideas are closely related 
to the different point of views on the origin of the pseu- 
dogap and its relation to the superconducting phase. 

For example, if one takes the point of view that the 
pseudogap is related to an order parameter which com- 
petes with superconductivity then the tendency to de- 
velop such an order reduces T c . Therefore, if in the het- 
erostructure this tendency is suppressed superconductiv- 
ity will prevail to higher temperatures. Some support for 
this point of view is provided by the experiment [l| in the 
following way. The largest enhancement occurs close to 
1/8 filling. At this filling, the single compound films have 
a dip in the T c vs. x curve. The dip is believed to be 
due to charge stripes which are favorable at this filling. 
In the heterostructures no such dip is observed. 

Another point of view of the pseudogap physics is that 



in the underdoped side of the superconducting dome, T c 
is restricted by phase fluctuations while pairing persists 
up to higher energy scale, possibly of the order of the 
pseudogap temperature, T*. Taking this point of view, it 
has been suggested by Yuli et al. [l[ and theoretically ex- 
plored by Berg et ai.[4| and by Goren and AltmanQ, that 
in the heterostructure, enhancement may occur if phase 
fluctuations are suppressed compared to the single com- 
pound film. To suppress the fluctuations the supcrfluid 
density should be increased, making vortices more costly 
in energy. Higher superfluid density can be achieved by 
higher density of states at the Fermi level which can be 
provided by a metallic layer. This view is supported by 
the fact that the enhancement only occurs on the under- 
doped side of the superconducting dome. 

Both of the directions above are interesting and deserve 
further investigation. In this work, however, we point at 
a different direction which does not necessarily discrimi- 
nate between the two scenarios above but may lead the 
way to a microscopic description of the phenomenon. We 
are guided by both experimental evidence and relevant 
numerical findings. We first assume that the cause for the 
T c enhancement resides within the copper-oxygen plane. 
This starting point should be contrasted with approaches 
that include a bilayer interface of the two materials and 
rely on strong hopping between the layers 0, [Bj]. Unfor- 
tunately, the inter-layer hopping strength in the cuprate 
materials may be too weak to explain such strong vari- 
ations in T c . A quantitative measure of the interlayer 
hopping strength has been seen, for example, in transport 
properties [fl 0] and optical conductivity]^. There is no 
evidence to suggest the heterostructures have stronger 
inter-plane coupling than that of the single compound 
film. 

What is then the difference between the single com- 
pound film and the heterostructure? One difference may 
be the electronic density. The amount of doping x in the 
bottom layer of the heterostructure may be different from 
that of the starting material due to charge migration be- 
tween the two layers. This leads to a 'self doping' effect 
which may change T c . This effect alone, however, cannot 
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account for the large T c enhancement in the experiment. 
Doping alone can not yield a transition temperature that 
is significantly higher than that of the single compound 
film at optimal doping. 

Another difference is structural; atoms move and 
bonds stretch/shrink to relieve strain resulting from the 
two layer mismatch. A recent experiment by Zhou et 
al.^ provides evidence that structural changes arc inti- 
mately related to the T c enhancement. In these exper- 
iments heterostructures of the parent compound (LCO) 
and overdoped metallic LSCO are fabricated. These het- 
erostructures also display superconductivity with an en- 
hanced T c relative to the single film. In order to de- 
tect where in the heterostructure superconductivity oc- 
curs, the experimental group introduce zinc impurities 
to the sample, layer by layer (loj. In-plane zinc impu- 
rities are known to suppress T c in the cuprates by a 
factor of about 2. In the heterostructure their effect 
was significant only when introduced to one specific layer 
which resides 1 layer above the interface, on the insulat- 
ing/superconducting LCO side. This leads to the conclu- 
sion that only one layer is responsible for the heterostruc- 
ture's high T c . In addition, the apical oxygen distance 
from each layer (d) was measured by x-ray scattering 
and was found to vary significantly in the heterostruc- 
ture, between d as 2.3A and d as 2.75A. In contrast, 
both the insulator (LCO) and metal (LSCO) have a bulk 
apical distance of d as 2. 4 A (as pointed out in Ref. Ill 
and measured in Ref.[H|). Therefore, the heterostructure 
achieves apical oxygen distances in doped LSCO, that are 
not reached in bulk samples. We believe that this lattice 
rearrangement is responsible for the T c enhancement and 
motivate this point of view using a microscopic model be- 
low. 

In order to discuss superconductivity in the cuprates 
without speculating on the pairing mechanism and the 
origin of the pseudogap we explore the extended Hub- 
bard model. This model contains both the hopping pa- 
rameters that are relevant to the cuprates and can be 
fit to the band structure as it is observed in experiment 
(e.g. ARPES 12|) and the strongly correlated nature of 
this electronic system. The Hamiltonian reads 

ij,a i (ij) 

where the indices i, j go over the lattice sites with angular 
brackets denoting nearest neighbors and a =t4 is the 
spin label. The number of electrons at site i is denoted 
by rii = riif + n ti = J2 a C L C ^- 

Numerical analysis of this model leads to d-wave 
superconductivity [13J. In order to see how this comes 
about one has to go beyond simple mean field and in- 
clude strong spin and charge fluctuations. In this model 
the onsite repulsion U promotes d-wave superconductiv- 
ity while the off-site repulsion V suppresses it. This 
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FIG. 1: The effective attractive off-site interaction V" e // as a 
function of the bare repulsive, off-site, interaction V in units 
of the on-site repulsion U. In the calculation the values t = 
488meV and U = 2.56eV were used. The (bare) repulsive 
V axis is in units of meV. The values taken for the bare 
repulsive V range from V — to its expected value for bulk 
LCO as found in Ref. GJ. 



has been seen, for example, in the fluctuation exchange 
(FLEX) approximation [13j. In this approximation spin 
and charge fluctuations are taken into account (to all 
orders in perturbation theory) in the renormalized in- 
teraction vertex function. This vertex together with the 
electronic Green's function is used in the Eliashberg[l3| 
theory to determine the strength of pairing. The analysis 
shows that when the off-site interaction V is increased 
d-wave pairing is suppressed. If V is larger than some 
critical value a charge density wave becomes more ener- 
getically favorable than superconductivity. 

The connection between the off-site interaction term 
and the apical oxygen distance from the Cu-0 plane has 
been recently explored from first principles by Yin and 
Ku[l5| and by Weber et a2.[l(|. The conclusion of these 
studies is that the distance of the apical oxygen from 
the layer is inversely related to the bare V. For mate- 
rials with larger apical distance the off-site repulsion is 
smaller. Other terms in the Hamiltonian, however, are 
rather insensitive to this distance. The result is a d- 
wave superconductor with larger amplitude. This is the 
main point of the paper - the apical oxygen distances in 
the heterostructure are significantly larger than in bulk 
LSCO, in particular at the superconducting layer. This 
causes a reduction of the off-site repulsion V and super- 
conductivity is therefore enhanced. 

In order to relate these findings to the experimental 
setup we investigate the Hamiltonian in Eq.[T]as a model 
for the highest T c layer close to the interface. We start 
with repulsive interactions both on-site (U) and off-site 
(V) and vary the strength of V to capture the effect of 
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FIG. 2: The mean field result for d-wave superconductivity 
order parameter (squares) and Neel order (diamonds) as a 
function of the ratio between the off-site interaction V and 
the on-site interaction U. The inset shows a zoom- in on the 
relevant interaction ratio for the superconducting order pa- 
rameter Aq. 



the apical oxygen distance. As a first step we use the 
FLEX approximation to include the Fermionic response 
in the effective interaction strength. The renormalized 
interactions U e ff and V e ff are presented in Fig.[TJ Note 
that the effective off-site interaction is attractive. Using 
renormalized interaction coefficients is crucial, since su- 
perconductivity requires pairing channel. 

An intuitive understanding of the mechanism by 
which off-site interaction is related to superconductiv- 
ity emerges from the FLEX approximation. We start 
with a repulsive on-site and off-site interaction in the ex- 
tended Hubbard model, as appropriate for this system. 
However, when considering both interactions together in 
the vertex function the total off-site channel is attrac- 
tive. This can be seen even at the level of the standard 
Hubbard model (without off-site interaction). The Hub- 
bard model close to half filling was studied by Scalapino 
in Ref. 



17j where it was found that the on-site repul- 



sive interaction U leads to an effective off-site attraction. 
This is the result of proximity to the antiferromagnetic 
phase at half filling which is due to double hopping pro- 
cesses of the order of t 2 /U. Away from half filling there 
is no long range Neel order but the susceptibility is still 
peaked around (tt, tt). The contribution of such a struc- 
ture to the interaction vertex is, again, strongest at this 
point in momentum space which leads to a large attrac- 
tion on nearest neighbor sites. This can be viewed as 
large Friedel oscillations. Electrons on nearest neighbor 
sites take advantage of these oscillations through pairing 
in the d-wave channel. 

In order to determine T c realistically one should use 
the renormalized vertex in a full Eliashberg calculation as 



was done by Onari et aZ.[13j] In the present work we cal- 
culate effective interaction parameters using the FLEX 
approximation, and turn to mean field theory in order 
to sketch the qualitative effect of the interaction on the 
order parameter. In the effective model the on-site in- 
teraction U e ff is repulsive while the off-site interaction 
V e ff is attractive. The bare repulsion, and the induced 
effective attraction combine to give an effective attractive 
off-site interaction and a renormalized on-site repulsion. 
When the apical oxygen distance is increased, the bare 
repulsion is decreased, and the overall effective interac- 
tion is more strongly attractive. We show the mean-field 
result for the above effective Hubbard model in Fig [21 In 
this figure, as the ratio between the magnitude of the at- 
tractive V e ff and the on-site interaction U e f / is increased 
d-wave superconductivity is enhanced and the Neel or- 
der (or spin density wave) is suppressed. This result is 
obtained by the usual mean field decomposition of the 
interaction terms in the superconducting and antiferro- 
magnetic channels. The order parameters are found by 
minimizing the mean-field free energy. 
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where RBZ stands for the sum over the reduced Brillouin 
zone, e k = — 2<(cos(fca;) + cos(fcj,)) — fi + Un is the band 
dispersion which includes a chemical potential shift Un, 
A k = Aox k with Xk = 2(cos(k x ) — cos(k y )) is the gap 
function, S is the Neel order parameter which folds the 
Brillouin zone with the wavevector Q = (tt, tt) and N is 
the number of sites. We have used an attractive potential 
to derive these equations and write |V e //| to emphasize 
this point. 
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where the third equation has only a trivial solution 
(ri = 0) at half filling. We must point out the numbers we 
find from our simple mean field treatment in Fig. [2] sug- 



gest the superconductor is significant only if ■ 
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> 0.2, 



while Fig. [T] shows a value no higher than 0.08 for this 
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ratio. However, in a full Eliashberg treatment even a 
small attractive off-site potential (like that of Fig. [IJ can 
cause the superconducting T c to be non-zero [lj|. We 
therefore conclude that our simplistic analysis gives the 
correct qualitative results, while leaving something to be 
desired in terms of quantitative results. 

Before our conclusions, let us briefly discuss the possi- 
bility that phase fluctuations are suppressing T c on the 
underdoped side. In the models above only the ampli- 
tude of the d-wave pairing order parameter is considered. 
However, there is a lot of experimental evidence that 
phase fluctuations are important in the pseudogap phase. 
Naturally, the question arises whether phase fluctuations 
are reducing T c from the temperature in which pairs are 
formed (which should be associated with a higher en- 
ergy scale, possibly as high as T*) or pairing and phase 
coherence appear at the same temperature (T c ) while 
the higher energy pseudogap is due to a competing order 
parameter. To shed more light on this issue it would 
be helpful to study the pseudogap regime of the het- 
erostructure and measure T* . We propose to perform the 
following experiments. Though experimentally challeng- 
ing, it may be possible to measure the density of states 
(DOS) in the heterostructure and determine its pseudo- 
gap temperature T* . If we adopt the point of view that 
T* is due to a competing order, unrelated to supercon- 
ductivity we expect this temperature to be lower in the 
heterostructure than in the single compound films. This 
will reflect the competition; when the competing order 
is suppressed by a larger apical oxygen distance, d-wave 
superconductivity is enhanced. If, on the other hand T* 
is the temperature at which pairing begins it should be 
higher in the heterostructure since T c and T* are related. 
Another important measurement is that of the Nernst ef- 
fect which has provided evidence for phase fluctuations 
in the past [3] together with a diamagnetism probe (lij. 
It is important to determine whether the phase fluctua- 
tions are enhanced or suppressed in the heterostructure 
relative to the single compound film. 

In addition, we would like to propose a third experi- 
ment which would directly address the role of the apical 
oxygen in the T c enhancement. Applying tensile stress 
along the c-axis will act to reduce the apical oxygen dis- 
tance. This, according to the scenario we present here 
should increase the bare off-site repulsion and reduce T c . 

To summarize, we provide a simple explanation to 
the enhancement of the superconducting transition tem- 
perature, T c , in LSCO heterostructures compared with 
heterogenous, single compound films, as seen in recent 
experiments P, Q. We have analyzed the extended Hub- 
bard model in the FLEX approximation and used a 
heuristic mean-field analysis to asses the effect of the 
apical oxygen distance on order formation. We con- 
clude that the off-site nearest neighbor interaction has 
a tremendous effect on T c in d-wave superconductors. In 
particular, when this interaction becomes more negative 



(in other words when it's bare repulsive component is re- 
duced) the transition temperature increases. This inter- 
action has been found to be inversely related to the apical 
oxygen distance from the copper oxygen planes 15. la ]. 
It has also been seen that the apical oxygen distance is 
larger in the heterostructure, compared to the single com- 
pound filmQ. We conclude that it is this structural dif- 
ference between the heterostructure and the single com- 
pound film that is responsible for the enhancement of 
T c . 
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